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Abstract 

Consider transportation of one distribution of mass onto another, chosen to optimize 
the total expected cost, where cost per unit mass transported from a; to ?/ is given by 
a smooth function c{x,y). If the source density f~^{x) is bounded away from zero and 
infinity in an open region U' C R", and the target density f~{y) is bounded away from 
zero and infinity on its support V C R", which is strongly c-convex with respect to J7', 
and the transportation cost c is non-negatively cross-curved, we deduce continuity and 
injectivity of the optimal map inside U' (so that the associated potential u belongs to 
C^{U')). This result provides a crucial step in the low/interior regularity setting: in a 
subsequent paper [15] , we use it to establish regularity of optimal maps with respect to 
the Riemannian distance squared on arbitrary products of spheres. The present paper 
also provides an argument required by Figalli and Loeper to conclude in two dimensions 
continuity of optimal maps under the weaker (in fact, necessary) hypothesis (A3w) 
[17] . In higher dimensions, if the densities are Holder continuous, our result permits 
continuous differentiability of the map inside U' (in fact, C^^^" regularity of the associated 
potential) to be deduced from the work of Liu, Trudinger and Wang [33j . 
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1 Introduction 

Given probability densities < G L^(R") with respect to Lebesgue measure .if" on R", 
and a cost function c : R" x R" i — > [0, +oo], Monge's transportation problem is to find a 
map G : R" i — > R" pushing d/^"*" = f'^d^^ forward to d^i" = f~d^^ which minimizes the 
expected transportation cost [38] 

inf _ [ c{x,G{x))dn{x), (1.1) 

where G^fi'^ = fi~ means /U^[y] = fi'^[G~^{Y)] for each Borel Y C R". 

In this context it is interesting to know when a map attaining this infimum exists; sufficient 
conditions for this were found by Gangbo [20] and by Levin [31], extending work of a number 
of authors described in [21] [46]. One may also ask when G will be smooth, in which case it 
must satisfy the prescribed Jacobian equation | det DG{x)\ = f~^{x)/f~{G{x)), which turns 
out to reduce to a degenerate elliptic partial differential equation of Monge- Ampere type for 
a scalar potential u satisfying Du{x) = —Dxc{x^G{x)). Sufficient conditions for this were 
discovered by Ma, Trudinger and Wang [ST] and Trudinger and Wang [43] [44], after results 
for the special case c{x,y) = \x — 2/P/2 had been worked out by Brenier [4], Delanoe |12| . 
Caffarelli [6] [5] [7] [8] [9], and Urbas [45], and for the cost c(x, y) = — log |x — 7/| and measures 
supported on the unit sphere by Wang [48]. 

If the ratio f^{x)/f~{y) — although bounded away from zero and infinity — is not con- 
tinuous, the map G will not generally be differentiable, though one may still hope for it to 
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be continuous. This question is not merely of academic interest, since discontinuities in / 
arise unavoidably in applications such as partial transport problems [lO] [3] [H] [14]. Such 
results were established for the classical cost c{x,y) = \x — by Caffarelli [5] [7] [8], for 

its restriction to the product of the boundaries of two strongly convex sets by Gangbo and 
McCann [22], and for more general costs satisfying the strong regularity hypothesis (A3) of 
Ma, Trudinger and Wang [37] — which excludes the cost c{x,y) = \x — yp/2 — by Loeper 
|34| : see also [27] [32] [44]. Under the weaker and degenerate hypothesis (A3w) of Trudinger 
and Wang [43], which includes the cost c{x, y) = \x — yp/2 (and whose necessity for regular- 
ity was shown by Loeper [31]), such a result remains absent from the literature; we aim to 
provide it below under a slight strengthening of their condition (still including the quadratic 
cost) which appeared in Kim and McCann [28] [29], called non-negative cross-curvature. (Re- 
lated but different families of strengthenings were investigated by Loeper and Villani [36] and 
Figalli and Rifford [IS].) Our main result is stated in Theorem 12.11 A number of interesting 
cost functions do satisfy non-negative cross-curvature hypothesis, and have applications in 
economics [16] and statistics [40]. Examples include the Euclidean distance between two con- 
vex graphs over two sufficiently convex sets in R*^ [37], the Riemannian distance squared on 
multiple products of round spheres (and their Riemannian submersion quotients, including 
products of complex projective spaces CP") [29], and the simple harmonic oscillator action 
|30| . In a sequel, we apply the techniques developed here to deduce regularity of optimal 
maps in the latter setting [15]. Moreover, Theorem 12.11 allows one to apply the higher in- 
terior regularity results established by Liu, Trudinger and Wang [33], ensuring in particular 
that the transport map is C°°-smooth if and /~ are. 

Most of the regularity results quoted above derive from one of two approaches. The 
continuity method, used by Delanoe, Urbas, Ma, Trudinger and Wang, is a time-honored 
technique for solving nonlinear equations. Here one perturbs a manifestly soluble problem 
(such as I detL'Go(x)| = f~^{x)/fQ{Go{x)) with fo = f~^, so that Gq{x) = x) to the problem 
of interest (|detDGi(x)| = f^{x)/ fi{Gi{x)), fi = /") along a family {ft}t designed to 
ensure the set of t £ [0, 1] for which it is soluble is both open and closed. Openness follows 
from linearization and non-degenerate ellipticity using an implicit function theorem. For the 
non-degenerate ellipticity and closedness, it is required to establish estimates on the size of 
derivatives of the solutions (assuming such solutions exist) which depend only on information 
known a priori about the data (c, /t). In this way one obtains smoothness of the solution 
y = Gi{x) from the same argument which shows Gi to exist. 

The alternative approach relies on first knowing existence and uniqueness of a Borel map 
which solves the problem in great generality, and then deducing continuity or smoothness 
by close examination of this map after imposing additional conditions on the data (c, /^). 
Although precursors can be traced back to Alexandrov [2], in the present context this method 
was largely developed and refined by Caffarelli [5] [7] [8], who used convexity of u crucially to 
localize the map G{x) = Du{x) and renormalize its behaviour near a point (x, G{x)) of inter- 
est in the borderline case c{x,y) = —{x,y). For non-borderline (A3) costs, simpler estimates 
suffice to deduce continuity of G, as in [22] [H] [S] [H]; in this case Loeper was actually 
able to deduce an explicit bound a = (4n — 1)~^ on the Holder exponent of G when n > 1, 
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which was recently improved to its sharp value a = (2n — by Liu [32] using a technique 
related to the one we develop below and discovered independently from us; both Loeper and 
Liu also obtained explicit exponents a = a{n,p) for G with p > n [34] or p > (n + 1)/2 
[32] and G L°°. Explicit bounds on the exponent are much worse in the classical case 

c{x,y) = —{x,y) [19], when such exponents do not even exist unless log G L°° [8] [47] . 

Below we extend the approach of Caffarelli to non-negatively cross-curved costs, a class 
which includes the classical quadratic cost. Our idea is to add a null Lagrangian term to the 
cost and exploit diffeomorphism (i.e. gauge) invariance to choose coordinates which depend 
on the point of interest that restore convexity of u{x); our strengthened hypothesis then 
permits us to exploit Caffarelli's approach more systematically than Liu was able to do [32]. 
However, we still need to overcome serious difficulties, such as getting an Alexandrov estimate 
for c-subdifferentials (see Section [T]) and dealing with the fact that the domain of the cost 
function (where it is smooth and satisfies appropriate cross-curvature conditions) may not 
be the whole of R". (This situation arises, for example, when optimal transportation occurs 
between domains in Riemannian manifolds for the distance squared cost or similar type.) 
The latter is accomplished using Theorem 15.11 where it is first established that optimal 
transport does not send interior points to boundary points, and vice versa, under the strong 
c-convexity hypothesis (B2u) described in the next section. ( For this result to hold, the 
cost needs not to satisfy the condition (A3w).) Without our strengthening of Trudinger and 
Wang's hypothesis [43] (i.e. with only (A3w)), we obtain the convexity of all level sets of 
u{x) in our chosen coordinates as Liu also did; this yields some hope of applying Caffarelli's 
method and the full body of techniques systematized in Gutierrez [23], but we have not 
been successful at overcoming the remaining difficulties in such generality. In two dimensions 
however, there is an alternate approach to establishing continuity of optimal maps which 
applies to this more general case; it was carried out by Figalli and Loeper [T^, but relies on 
Theorem 15.11 first proved below. 



2 Main result 

Let us begin by formulating the relevant hypothesis on the cost function c{x, y) in a slightly 
different format than Ma, Trudinger and Wang [37]. For each {x^y) ^ U x V assume: 
(BO) ?7 C R" and y C R" are open and bounded and c G ^^(17 x F); 

\ n • , • , N X e U I — > —Dycix, w) 1 , . , ,1 . 

(Bl) bi-twist -TT )~ { r are difteomorphisms onto their ranges; 

^ ^ ^ ^ y G y I — > -D^c{x,y) J ^ ^ ' 

(B2) (bi-convex) ^/ '._ n^'^i^^f)\ \ are convex subsets of R" ; 
(B3) (non-negative cross-curvature) 



& 

cross(^(o),j^(o))[x'(0),y'(0)] := 



c{x{s),y{t))>0 (2.1) 

(s,t)={0,0) 
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for every curve t € [—1,1] i — > {Dyc{x{t),y{0)), Dxc{x{0),y{t))) G R^" which is an affinely 
parameterized line segment. 

If the convex domains Uy and 14 in (B2) are all strongly convex, we say (B2u) holds. 
Here a convex set U C is said to be strongly convex if there exists a radius R < +00 
(depending only on [/,) such that each boundary point x E dU can be touched from outside 
by a sphere of radius R enclosing U; i.e. U C Br{x — Rnu{x)) where nu{x) is an outer 
unit normal to a hyperplane supporting U at x. When U is smooth, this means all principal 
curvatures of its boundary are bounded below by 1/R. Hereafter U denotes the closure of 
[/, int U denotes its interior, diamf/ its diameter, and for any measure /i^ > on U , we use 
the term support and the notation spt //"'" C Z7 to refer to the smallest closed set carrying the 
full mass of ^l'^ . 

Condition (B3) is the above-mentioned strengthening of Trudinger and Wang's criterion 
(A3w) guaranteeing smoothness of optimal maps in the Monge transportation problem (11.11) : 
unlike us, they require (12. ip only if, in addition |43| . 

92 



dsdt 



c(x(s),y(t)) =0. (2.2) 

(s,t)={0,0) 



Necessity of Trudinger and Wang's condition for continuity was shown by Loeper [3l], who 
noted its covariance (as did [28] [H]) and some relations to curvature. Their condition 
relaxes the hypothesis (A3) proposed earlier with Ma [3^, which required strict positivity 
of (|2.ip when (12. 2p holds. The strengthening considered here was first studied in a different 
but equivalent form by Kim and McCann [28], where both the original and the modified 
conditions were shown to correspond to pseudo-Riemannian sectional curvature conditions 
induced by the cost c on ?7 x y, highlighting their invariance under reparametrization of 
either U 01 V hy diffeomorphism; see [28l Lemma 4.5]. The convexity of Uy required in 
(B2) is called c-convexity of U with respect to y by Ma, Trudinger and Wang (or strong 
c-convexity if (B2u) holds); they call curves x{s) £ U , for which s £ [0, 1] 1 — > —Dyc{x{s), y) 
is a line segment, c-segments with respect to y. Similarly, V is said to be strongly c*-convex 
with respect to x — or with respect to U when it holds for all x £ U — and the curve 
y{t) from (|2.1|) is said to be a c*-segment with respect to x. Such curves correspond to 
geodesies {x{t),y) and {x,y{t)) in the geometry of Kim and McCann. Here and throughout, 
line segments are always presumed to be affinely parameterized. 
We are now in a position to summarize our main result: 

Theorem 2.1 (Interior continuity and injectivity of optimal maps). Let c S C^([/ x V) 
satisfy (B0)-(B3) and (B2u). Fix probability densities /+ G L^{U) and f~ £ L^i^) with 
(/^//") G L'^^iU X V) and set dfi^ := // the ratio {f / f+) G L°°(C/' ^V) for 

some open set U' C U, then the minimum ^-1]) is attained by a map G : U 1 — > V whose 
restriction to U' is continuous and one-to-one. 

Proof. As recalled below in Section[3](or see e.g. [46]) it is well-known by Kantorovich duality 
that the optimal joint measure 7 G T{fj.~^, from (|3.ip vanishes outside the c-subdifferential 
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(|3.3p of a potential u = u'^*'^ satisfying the c-convexity hypothesis (|3.2I) . and that the map 
G : U i — > V which we seek is uniquely recovered from this potential using the diffeomorphism 
(Bl) to solve (13. Sp . Thus the continuity claimed in Theorem 12. II is equivalent to u S C^{U'). 

Since fi^ do not charge the boundaries of U (or of V), Lemma ISTT e) shows the c-Monge- 
Ampere measure defined in (|3.6I) has density satisfying \d'^u\ < ||/'''//~||L°°(t/xy) on U and 

\\r/f^\\lL(U'xV) ^ l^'^l ^ Wf^/rhooiu'xv) on U'. Thus u e C\U') according to 
Theorem 19.21 Injectivity of G follows from Theorem 19.11 and the fact that the graph of G is 
contained in the set S'^u C 17 X F of ([33]) • □ 

Note that in case G Cc{U) is continuous and compactly supported, choosing U' = 
U'^ = {f~^ > e} for all e > 0, yields continuity and injectivity of the optimal map y = G{x) 
throughout U'q. 

Theorem 12.11 provides a necessary prerequisite for the higher interior regularity results 
established by Liu, Trudinger and Wang in [33] — a prerequisite which one would prefer 
to have under the weaker hypotheses (B0)-(B2) and (A3w). Note that these interior 
regularity results can be applied to manifolds, after getting suitable stay-away-from-the-cut- 
locus results: this is accomplished for multiple products of round spheres in [15], to yield the 
first regularity result that we know for optimal maps on Riemannian manifolds which are not 
flat, yet have some vanishing sectional curvatures. 



3 Bacltground, notation, and preliminaries 

Kantorovich discerned [25] [26] that Monge's problem (II. ip could be attacked by studying 
the linear programming problem 

min / c{x,y)d'y{x,y). (3.1) 

'ye^{^i+,^i )JuxV 

Here T{fi'^,^^) consists of the joint probability measures on U x V C R" x R" having //^ 
for marginals. According to the duality theorem from linear programming, the optimizing 
measures 7 vanish outside the zero set of u{x) + v{y) + c(x, y) > for some pair of functions 
(u, v) = {v^, ) satisfying 

v''{x) ■.= s\yp-c{x,y) -v{y), (y) := sup -c(x, y) - n(x); (3.2) 

y&V xeU 

these arise as optimizers of the dual program. This zero set is called the c-subdifferential of 
u, and denoted by 

8% = {{x,y) £U xV \ u{x) + u"* (y) + c{x, y) = 0}; (3.3) 

we also write d'^u{x) := {y \ {x,y) G d'^u}, and d'^*u'^*{y) := {x \ (x,y) G d'^u}, and 
d'^u{X) := Ux&xd'^u{x) for X C R". Formula (|3.2p defines a generalized Legendre-Fenchel 
transform called the c-transform; any function satisfying u = u'^ := {u^ Y is said to be 
c-convex, which reduces to ordinary convexity in the case of the cost c{x,y) = —{x,y). In 
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that case d'^u reduces to the ordinary subdifferential du of the convex function u, but more 
generally we define 

du := {{x,p) £ U X R" I u{x) > u{x) + {p, x — x) + o{\x — x\) as x — > x}, (3.4) 

du{x) := {p I {x,p) G du}, and du{X) := Llx£xdu{x). Assuming c G C^[U x V) (which 
is the case if (BO) holds), any c-convex function u = u'^ will be semi-convex, meaning its 
Hessian admits a bound from below D^u > — ||c||c2 in the distributional sense; equivalently, 
u{x) + ||c||c2 |xp/2 is convex on each ball in U [21| . In particular, u will be twice-differentiable 
^^-a.e. on U in the sense of Alexandrov. 

As in [20] [31] [37], hypothesis (Bl) shows the map G : dom Du i — > V is uniquely defined 
on the set dom Du C U of differentiability for u by 

Dxc{x,G{3:)) = -Du{x). (3.5) 

The graph of G, so-defined, lies in d'^u. The task at hand is to show continuity and injectivity 
of G — the former being equivalent to u G C^{U) — by studying the relation d'^u C U x V. 
To this end, we define a Borel measure \d^u\ on R" associated to u by 

la^nKAT) :=^"(a"n(X)) (3.6) 

for each X C R"; it will be called the c-Monge- Ampere measure of u. (Similarly, we define 
\du\.) We use the notation \d'^u\ > A on ?7' as a shorthand to indicate \d'^u\{X) > A^"(X) 
for each X C U'; similarly, \d^u\ < A indicates \d^u\{X) < A^"(X). As the next lemma 
shows, uniform bounds above and below on the marginal densities of a probability measure 
7 vanishing outside d'^u imply similar bounds on \d'^u\. 

Lemma 3.1 (Properties of c-Monge- Ampere measures). Let c satisfy (BO)-(Bl), while u 
and Uk denote c-convex functions for each /c G N. Fix x £ X and constants A, A > 0. 

(a) Then d^u{U) C V and \d^u\ is a Borel measure of total mass ^^[V) on U. 

(b) If Uk Uoo uniformly, then u^q is c-convex and \d^Uk\ l^^^^^ool weakly-* in the duality 
against continuous functions on U x V . 

(c) Ifuk{x) = for all k, then the functions u^ converge uniformly if and only if the measures 
\d'^Uk\ converge weakly-*. 

(d) If \d''u\ <k onU, then |a^*n^* \>l/K onV. _ _ 

(e) If a probability measure 7 > vanishes outside d'^u <ZU xV , and has marginal densities 
f^, then f + > X on U' CU and f~<AonV imply {d^ul > A/A on U' , whereas /+ < A 
on U' and f~>XonV imply \d'^u\ < A/A on U' . 

Proof, (a) The fact d'^u{U) C F is an immediate consequence of definition (13. 3|) . Since 
c £ C^{U xV), the c-transform v = u'^ : V 1 — > R defined by (13. 2p can be extended to a 
Lipschitz function on a neighbourhood of V, hence Rademacher's theorem asserts dom Dv is a 
set of full Lebesgue measure in V. Use (Bl) to define the unique solution F : domDv 1 — > U 
to 

Dyc{F{y),y) = -Dv{y). 
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As in [20] [M], the vanishing of u{x) + v{y) + c{x,y) > imphes d'^*v{y) = {F{y)}, at least 
for all points y S domDv where V has Lebesgue density greater than one half. For Borel 
X C R", this shows d'^uiX) differs from the Borel set F~^{X) n F by a negligible subset 
of V, whence \d'^u\ = F^[^"-[yj so claim (a) of the lemma is established. 

(b) Let \\uk — lioollioo^jj) — > 0. It is not hard to deduce c-convexity of Uoo, as in e.g. |16| . 
Define Vk = u'j^ and on domDvk C V as above, so that |f?'^tifc| = i^fc#(^'^[y). Moreover, 
Vk — > Voo in L°°{V), where Voo is the c*-dual to Uoo- The uniform semiconvexity of (i.e. 
convexity of Vk{y) + |||c||c2|?/P) ensures pointwise convergence of Dvk Dvoo ^"-a.e. on 
V. From Dyc{Fk{y),y) = —Dvk{y) we deduce F^ F^o =Sf"-a.e. on V. This is enough 
to conclude |9'^itfc| ^ \d'^Uk\, by testing the convergence against continuous functions and 
applying Lebesgue's dominated convergence theorem. 

(c) To prove the converse, suppose is a sequence of c-convex functions which vanish 
at X and ^ fioo weakly-*. Since the have Lipschitz constants dominated by ||c||(7i 
and U is compact, any subsequence of the Uk admits a convergent further subsequence by the 
Ascoli-Arzela Theorem. A priori, the limit Uoc might depend on the subsequences, but (b) 
guarantees |9'^Uoo| = /^oo, after which pH Proposition 4.1] identifies Uoo uniquely in terms of 
A^"*" = /^oo and n~ = Jf"[y, up to an additive constant; this arbitrary additive constant is 
fixed by the condition Uooix) = 0. Thus the whole sequence Uk converges uniformly. 

(e) Now assume a finite measure 7 > vanishes outside d^u and has marginal densities 
Then the second marginal dfi~ := f~d^^ of 7 is absolutely continuous with respect 
to Lebesgue and 7 vanishes outside the graph F : V 1 — > U , whence 7 = (F x id)^^" 
by e.g. [U Lemma 2.1]. (Here id denotes the identity map, restricted to the domain dom-Df 
of definition of F.) Recalling that \d'^u\ = F^[^^[yj (see the proof of (a) above), for any 
Borel X C U' we have 

A|9"n|(X) = A^"(F-^(X)) < / f-{y)d^''{y) = I f+{x)d^''{x) < A^"(A:) 

Jf-^(x) Jx 

whenever A < /" and < A. We can also reverse the last four inequalities and interchange 
A with A to establish claim (e) of the lemma. 

(d) The last point remaining follows from (e) by taking 7 = (F x ic?)^^". Indeed an 
upper bound A on \d'^u\ = F^^'^ throughout U and lower bound 1 on Jf" translate into a 
lower bound 1/A on u'^ |, since the reflection 7* defined by ^*{Y x X) := -flX x Y) for 
each X X Y C U X V vanishes outside d'^ and has second marginal absolutely continuous 
with respect to Lebesgue by the hypothesis \d^u\ < A. □ 

Remark 3.2 (Monge- Ampere type equation). Differentiating (|3.5p formally with respect to 
X and recalling | detDG(i;)| = /^(x)//~(G(x)) yields the Monge- Ampere type equation 

|detD2^c(x,G(x))| f~{G{x)) 

on U, where G{x) is given as a function of x and Du{x) by (13.51) . Degenerate ellipticity 
follows from the fact that y = G{x) produces equahty in u{x) + u'^ (y) + c{x, y) > 0. A 
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condition under which c-convex weak-* solutions are known to exist is given by 

/ /+(x)d^"(x)= [_r{y)d^^{y). 
Ju Jv 

The boundary condition d^u{U) C V which then guarantees Du to be uniquely determined 
/"•"-a-e. is built into our definition of c-convexity. In fact, [Ml Proposition 4.1] shows u to 
be uniquely determined up to additive constant if either /"'' > or > ^"-a.e. on its 
connected domain, U or V . 

A key result we shall exploit several times is a maximum principle first deduced from 
Trudinger and Wang's work [43] by Loeper; see [34l Theorem 3.2]. A simple and direct proof, 
and also an extension can be found in [2H Theorem 4.10], where the principle was also called 
'double-mountain above sliding-mountain' (DASM). Other proofs and extensions appear in 
[44] [42] [46] [36] [B]: 

Theorem 3.3 (Loeper 's maximum principle 'DASM'). Assume (B0)-(B2) and (A3w) 

and fix x^x £U. If t £ [0,1] i — > —Dxc{x,y{t)) is a line segment then f{t) := —c{x,y{t)) + 
c{x,y{t)) < max{/(0),/(l)} for all t e [0, 1]. 

It is through this theorem and the next that hypothesis ( A3w) and the non-negative cross- 
curvature hypothesis (B3) enter crucially. Among the many corollaries Loeper deduced from 
this result, we shall need two. Proved in [HU Theorem 3.1 and Proposition 4.4] (alternately 
[281 Theorem 3.1] and [27l A. 10]), they include the c-convexity of the so-called contact set 
(meaning the c*-subdifferential at a point), and a local to global principle. 

Corollary 3.4. Assume (B0)-(B2) and (A3w) and fix {x,y) £ U xV. If u is c-convex 
then d'^u{x) is c* -convex with respect to x £ U , i.e. —Dxc{x,d^u{x)) forms a convex subset 
of T^U . Furthermore, any local minimum of the map x £ U i — > u{x) -\- c{x,y) is a global 
minimum. 

As shown in [29j Corollary 2.11], the strengthening (B3) of hypothesis (A3w) improves 
the conclusion of Loeper's maximum principle. This improvement asserts that the altitude 
f{t,x) at each point of the evolving landscape then accelerates as a function of t G [0, 1]: 

Theorem 3.5 (Time-convex DASM). Assume (B0)-(B3) and fix x,x £U. If t £ 

[0, 1] I — > —Dxc{x, y{t)) is a line segment then t £ [0, 1] i — > f{t) := — c(x, y{t)) + c(x, y{t)) is 
convex. 

Remark 3.6. Since all assumptions (B0)-(B3) and (A3w) on the cost are symmetric in x 
and y, all the results above still hold when exchanging x with y. 

4 Cost-exponential coordinates, null Lagrangians, and afRne 
renormalization 

In this section, we set up the notation for the rest of the paper. Recall that c £ C^{U x V) 
is a non-negatively cross-curved cost function satisfying (B1)-(B3) on a pair of bounded 
domains U and V which are strongly c-convex with respect to each other (B2u). 
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Fix A, A > and an open domain U C U, and let u be a c-convex solution of the 
c-Monge- Ampere equation 

r X^"" < la'^nl < in U^CU, 

{ Id^ul < A^" in U. ^ ' 

We sometimes abbreviate (14. Ij) by writing \d'^u\ E [A, 1/A]. In the following sections, we will 
prove interior differentiability of u on ^7^, that is u G C^(C/^); see Theorem [921 

Throughout Dy will denote the derivative with respect to the variable y, and iterated 
subscripts as in D^y denote iterated derivatives. We also use 

Pt = PtiUxV):= \\iDlyC)^'U^^uxV) (4.2) 
7± = jt{U xV):=\\det{Dlyc)^^\\L^(^u^v) (4-3) 

to denote the bi-Lipschitz constants of the coordinate changes (14. 4p and the Jacobian 
bounds 7^ for the same transformation. Notice Jclc ^ 1 ^^r any cost satisfying (Bl), 
and equality holds whenever the cost function c{x, y) is quadratic. So the parameter ^^Ic 
crudely quantifies the departure from the quadratic case. The inequality j3'^ I3~ > 1 is much 
more rigid, equality implying D'^yC{x,y) is the identity matrix, and not merely constant. 

4.1 Choosing coordinates which convexify c-convex functions 

In the current subsection, we introduce an important transformation (mixing dependent and 
independent variables) for the cost c{x,y) and potential u{x), which plays a crucial role 
in the subsequent analysis. This change of variables and its most relevant properties are 
encapsulated in the following definition and theorem. In the sequel, whenever we use the 
expression c{q, •) or u{q) we refer to the modified cost function and convex potential defined 
here, unless otherwise stated. Since properties (B0)-(B3), (A3w) and (B2u) were shown 
to be tensorial in nature (i.e. coordinate independent) in [28] [34], the modified cost c inherits 
these properties from the original cost c with one exception: (14. 5|) defines a diffeomorphism 
q & Uy I — > x{q) € C/, so the cost c G C^{Uy x V) may not be smooth. However, its 
definition reveals that we may still differentiate c four times as long as no more than three 
of the four derivatives fall on the variable g, and it leads to the same geometrical structure 
(pseudo-Riemannian curvatures, including (12. ip ) as the original cost c since the metric tensor 
and symplectic form defined in [28] involve only mixed derivatives D\yC, and therefore remain 
functions of the coordinates (g, y) G Uy x V . 

Definition 4.1 (Cost-exponential coordinates and apparent properties). Given c G C^{U x 
V) strongly twisted (BO)-(Bl), we refer to the coordinates {q,p) Uy x Vx defined by 

q = q{x) = -Dyc{x,y), p = p{y) = -Dxc{x,y), (4.4) 

as the cost exponential coordinates from y & V and x & U respectively. We denote the 
inverse diffeomorphisms by x : Uy C T^V i — > U and y : Vx C T^U i — > V ; they satisfy 

q = -Dyc{x{q),y), p = -Dxc{x, y{p)). (4.5) 
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The cost c{q,y) = c{x{q),y) — c{x{q),y) is called the modified cost at y. A subset of U 
or function thereon is said to appear from y to have property A, if it has property A when 
expressed in the coordinates q£Uy. 

Remark 4.2. Identifying the cotangent vector 0©^ with the tangent vector Q®0 to C/ xV 
using the pseudo-metric of Kim and McCann [28] shows x{q) to be the projection to U of 
the pseudo-Riemannian exponential map exp^^ Q © 0; similarly y{p) is the projection to V 
of exp(j © P. Also, x{q) =: c*-expy q and y{p) =: c-exp^p in the notation of Loeper |34| . 

Our first contribution is the following theorem. For a non-negatively cross-curved cost 
(B3), it shows that any c-convex potential appears convex from y £ V. Even if the cost 
function is weakly regular (A3w), the level sets of the c-convex potential appear convex 
from y, as was discovered independently from us by Liu [32], and exploited by Liu with 
Trudinger and Wang [33]. Note that although the difference between the cost c{x,y) and the 
modified cost c{q, y) depends on y, they differ by a null Lagrangian c(x, y) which — being 
independent of y E F — does not affect the question of which maps G attain the infimum 
(jl.ip . Having a function with convex level sets is a useful starting point, since it enables 
us to apply Caffarelli's affine renormalization of convex sets approach and a full range of 
techniques from Gutierrez [23] to address the regularity of c-convex potentials. 

Theorem 4.3 (Modified c-convex functions appear convex). Let c £ [U x V) satisfying 
(B0)-(B2) be weakly regular (A3w). If u = u'^ ^ is c-convex on U, then u{q) = u{x{q)) + 
c{x{q),y) has convex level sets, as a function of the cost exponential coordinates q S Uy from 
y G V. If, in addition, c is non-negatively cross-curved (B3) then u is convex on Uy. In 
either case u is minimized at qo if y £ d'^u{x{qo)). Furthermore, u is c-convex with respect to 
the modified cost c{q,y) := c{x{q),y) — c{x{q),y) on Uy x V , and d^u{q) = d^u{x{q)) for all 

q e Uy. 

Proof. The final sentences of the theorem are elementary: c-convexity u = u'^ asserts 

u{x) = sup —c{x,y) — u'^ (y) and u'^ (y) = sup —c{x{q),y) — u{x{q)) = (y) 
yeV geUy 

from (|3.2I) . hence 

u{q) = sup -c{x{q),y) +c{x{q),y) -u^'iy) 
ydV 

= sup -c{q,y) -u^'iy), 
y&v 

and d^u{q) = d^u{x{q)) since all three suprema above are attained at the same y £ V. Taking 
y = y reduces the inequality u{q) + vP {y) -\- c{q, y) > to u{q) > —u'^ {y) , with equality 
precisely if y G d^u{q). It remains to address the convexity claims. 

Since the supremum u{q) of a family of convex functions is again convex, it suffices to 
establish the convexity of g G Uy \ — > —c{q, y) for each y £ V under hypothesis (B3). For a 
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similar reason, it suffices to estabfish the level-set convexity of the same family of functions 
under hypothesis (A3w). 

First assume (A3w). Since 

Dyc{q,y) = Dyc{x{q),y) := -q (4.6) 

we see that c-segments in Uy with respect to y coincide with ordinary line segments. Let 
g(s) = {l — s)qQ+sqi be any line segment in the convex set Uy. Define /(s, y) := —c{q{s), y) = 
—c{x{q{s)),y) + c{x{q{s)),y). Loeper's maximum principal (Theorem 13.31 above, see also 
Remark I3.6P asserts f{s,y) < max{/(0, y), /(I, y)}, which implies convexity of each set 
{q £ Uy \ —c{q,y) < const}. Under hypothesis (B3), Theorem 13.51 goes on to assert 
convexity of s G [0, 1] i — > f{s,y) as desired. □ 

The effect of this change of gauge on Jacobian inequalities is summarized in a corollary: 

Corollary 4.4 (Transformed c-Monge- Ampere inequalities). Using the hypotheses and no- 
tation of Theorem\4^ if {d^'ul G [A, A] C [0, oo] on U' C U, then \d^ii\ £ [A/7+,A7~] on 
U'- = -Dyc{U',y), where = -ft{U' x V) and = (3t{U' x V) are defined in M -M . 
Furthermore, 7^ := 7^(C/^ x V) < 7+7" and (3f := l3f{U'y y.V)< PtPc- 

Proof. Prom the Jacobian bounds \ det Dxq{x)\ G [l/7cr>7^] on U' , we find ^"(X)/7^ < 
^"(g(^)) < 7+^"(X) for each X C U'. On the other hand, Theorem HJ asserts d^u{q{X)) = 
d'^u{X), so the claim \d'^u\ G [A/7^,A7~] follows from the hypothesis \d'^u\ G [A, A], by def- 
inition (13. 6|) and the fact that q : U — > Uy from (|4.4p is a diffeomorphism; see (Bl). The 
bounds 7^ < J^Jc ^iid (3f < P^P^ follow from D'^yc{q,y) = D'^yC{x{q),y)Dqx{q) and 
Dqx{q) = -Dlyc{x{q),yy^. □ 

4.2 AfSne renormalization 

The renormalization of a function u by an affine transformation L : R" — > R" will be useful 
in Section [7] to prove our Alexandrov type estimates. Let us therefore record the following 
observations. Define 

u*{q) = \detL\-^/''u{Lq). (4.7) 

Here detL denotes the Jacobian determinant of L, i.e. the determinant of the linear part of 
L. 

Lemma 4.5 (Affine invariance of c-Monge- Ampere measure). Assuming (BO)-(Bl), given 
a c- convex function u : Uy 1 — > R and affine bijection L : R" i — > IU\ define the renormalized 
potential u* by (|4.7p and renormalized cost 

c,{q, y) = I det Lr2/«c(Lg, L*y) (4.8) 
using the adjoint L* to the linear part of L. Then, for all Borel Z C Uy, 

\du*\{L-'^Z) = \deiL\-'^\du\{Z), (4.9) 

\d^*u*\{L-^Z) = IdetLl-V^-^K^)- (4.10) 



12 



Proof. Prom (fO) we see p € du{q) if and only if | det G du*{L~^q), thus ([49]) 

follows from du*{L~^Z) = \ det L\-^/''L*{du{Z)) . Similarly, since ^ yields {u*Y*{y) = 
I detL|~^/"'M'^*(L*y), we see y G d^u{q) is equivalent to | det L|~^/"L*y G d^*u*{L~^q) from 
(lOl (and Theorem [43]), whence d^*u*{L-^Z) = \ detL|-2/"L*(a^n*(Z)) to establish glO]). 

□ 

As a corollary to this lemma, we recover the affine invariance not only of the Monge- 
Ampere equation satisfied by u{q) — but also of the c-Monge- Ampere equation it satisfies 
— under coordinate changes on V (which induce linear transformations L on T~V and L* 
on TyV): for q e Uy, 

d\du*\,-^_, , d\du\ , , , d|(9'^*-u*| , dld'^ul , , 

5 Strongly c-convex interiors and boundaries not mixed by d^u 

The subsequent sections of this paper are largely devoted to ruling out exposed points in 
Uy of sets on which ordinary convexity of the c-convex potential from Theorem 14.31 fails to 
be strict. This current section rules out exposed points on the boundary of Uy. We do this 
by proving an important topological property of the (multi-valued) mapping d^u C U x V. 
Namely, we show that the subdifferential d^u maps interior points of spt|5'^u| C U only to 
interior points of V, under hypothesis (|4.ip . and conversely that d'^u maps boundary points 
of U only to boundary points of V. This theorem may be of independent interest, and was 
required by Figalli and Loeper to conclude their continuity result concerning maps of the 
plane which optimize (A3w) costs [T7] . 

This section does not use the cross-curvature condition (B3) (nor A3w) on the cost 
function c € C^{U x V), but relies crucially on the strong c-convexity (B2u) of its domains 
U and V (but importantly, not on sptl^^^ul). No analog for Theorem 15.11 was needed by 
Caffarelli to establish C^'^ regularity of convex potentials u{x) whose gradients optimize the 
classical cost c(x, y) = —{x,y) [8], since in that case he was able to take advantage of the fact 
that the cost function is smooth on the whole of R" to chase potentially singular behaviour 
to infinity. (One general approach to showing regularity of solutions for degenerate elliptic 
partial differential equations is to exploit the threshold-hyperbolic nature of the solution 
to try to follow either its singularities or its degeneracies to the boundary, where they can 
hopefully be shown to be in contradiction with boundary conditions; the degenerate nature 
of the ellipticity precludes the possibility of purely local regularizing effects.) 

Theorem 5.1 (Strongly c-convex interiors and boundaries not mixed by d^u). Let c satisfy 
(BO)-(Bl) and u = he a c-convex function (which implies d^u{U) = V ), and A > 0. 

(a) If \d^u\ > \ on X C U and V is strongly c* -convex with respect to X, then interior 
points of X cannot he mapped hy d^u to houndary points ofV: i.e. {X x dV) D d^u C 
{dX X dV). 



13 



(b) If \d'^u\ < A on U , and U is strongly c-convex with respect to V , then boundary points 
of U cannot be mapped by d^u into interior points of V : i.e. dU 'xV is disjoint from 
d^u. 

Proof. Note that when X is open the conclusion of (a) implies d^u is disjoint from X x dV . 
We therefore remark that it suffices to prove (a), since (b) follows from (a) exchanging the 
role X and y and observing that \d^u\ < A implies \d'^ u'^ \ > 1/ A as in Lemma [XlT d). 

Let us prove (a). Fix any point x in the interior of X, and y G d^u{x). Assume by 
contradiction that y E dV. At {x, y) we use (BO)-(Bl) to define cost-exponential coordinates 
I — ' {x{q),y{p)) by 

p= -D^c{x,y{p)) + D^c{x,y) ^T^{U) 

1= DlyC{x,y)-'^{Dyc{x{q),y) - Dyc{x,y)) <^Ti{U) 

and define a modified cost and potential by subtracting null Lagrangian terms: 

c{q,p) := c{x{q),y{p)) - c{x{p),y) - c{x,y{p)) 
u{q) ■= u{x{q)) + c{x{q),y). 

Similarly to Corollary 14.41 \d'^u\ > A := A/(7^"7~), where 7^ denote the Jacobian bounds 
(|4.3p for the coordinate change. Note = (a;(0),y(0)) corresponds to {p,q) = (0,0). 

Since c-segments with respect to y correspond to line segments in Uy := —Dyc{U,y) we see 
Dpc{q, 0) depends linearly on g, whence DggpC{q, 0) = 0; similarly c*-segments with respect 
to X become line segments in the p variables, Dqc{0,p) depends linearly on p, Dppgc{0,p) = 0, 
and the extra factor D'^yC{x,y)~^ in our definition of x{q) makes —Dp C{0,0) the identity 
matrix (whence q = —Dpc{0, q) and p = —Dqc{p, 0) for all q in Uy = x {U) and p in Vx := 
y~^{V)). Although the change of variables {q,p) 1 — > {x{p), y(?)) is only a diffeomorphism, 
we can still take four derivatives of the modified cost provided at least one of the four 
derivatives is with respect to q and another is with respect to p. We denote Xy := x"^(X) 
and choose orthogonal coordinates on U which make —in the outer unit normal to Vx C TfU 
at p = 0. Note that Vx is strongly convex by hypothesis (a). 

In these variables, consider a small cone of height e and angle 9 around the —in axis: 



{<? G 1 


- e - — 








kl 







Observe that, if 6,e are small enough, then Egj. C Xy, and its measure is of order e^6"' . 
Consider now a slight enlargement 

Ee,Cos ■■= {p = {P,Pn) G R." I Pn < e\p\ + Coe\p\'}, 

of the polar dual cone, where e will be chosen sufficiently small depending on the large 
parameter Co forced on us later. 

The strong convexity ensures Vx is contained in a ball BR{Rin) of some radius R > 1 
contained in the half-space Pn > with boundary sphere passing through the origin. As long 
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Figure 1: If d'^u sends an interior point onto a boundary point, then by c-monotonicity of d'^u the 
small cone Eg^^ has to be sent onto Eg n V^. Since for e > small but fixed ^"{Eg^^) ^ 
while J!f"{E'g(j^^ n Vx) < 0"''^^ (by the uniform convexity of Vj), we get a contradiction as 



as CqE < (6i?) ^ we claim Eg intersects this ball — a fortiori Vx — in a set whose volume 
tends to zero like 9^~^^ as 9 ^ 0. Indeed, from the inequality 

Pn<eVW+^+]:\P\ + ^Pn 

6 

satisfied by any {P,Pn) G £^e,Coe ^ Bn^Ren) we deduce pi < + 96'2)/(2 - 99^), i.e. 

Pn < \P\ 9 is small enough. Combined with the further inequalities 

^ < Pn < 0^\P\^+pI + Coe\P\^ + Coepl 

(the first inequality follows by the strong convexity of Vx), this yields \P\ < ^9\f2 and 
Pn < 0(9'^) as 6* ^ 0. Thus ^""{E'g^^^ n Vx) < C(9"+i for a dimension dependent constant C, 
provided CqE < (6i?)-^ 

The contradiction now will come from the fact that, thanks to the c-cyclical monotonicity 
of d'^u, if we first choose Cq big and then we take e sufficiently small, the image of all q G E'g^e 

by d'^u has to be contained in Eg ^^^^ for 9 small enough. Since d^u(^Xy^ C Vx this will imply 

~ ^^''{Ee^e) < \d'u\{Ee,e) < i^"(^x n E'g^Cos) < C9''+\ 

which gives a contradiction as ^ ^ 0, for e > small but fixed. 

Thus all we need to prove is that, if Co is big enough, then d'^u{E9^^) C Eg for any e 
sufficiently small. Let q G Eq^^ and p £ d'^u{q). Combining 

[ ds [ dt Dl^cisq, tp) [q,p\= c{q, p) + c(0, 0) - c(g, 0) - c(0, p) < 
JO JO 
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(where the last inequality is a consequence of c-monotonicity of d^u; see for instance [461 
Definitions 5.1 and 5.7]) with 

Dl c{sq, tp) =D^ c{0, tp) + I ds'Dl c{s'q, tp) [q] 

Jo 

=D^gpC{0,0)+ [ dt'DlpC{0,t'p)\p] 
Jo 

+ / ds'D^gpC{s'q,0)[q]+ / ds' / dt' D^gppC{s'q,t'p)[q,p] 
Jo Jo Jo 



yields 



hP) < 



/ ds dt ds' dt'D^gppC{s'q,t'p)[q,q,p,p] 
Jo Jo Jo Jo 



2|^|2 



< Co\q\'\p 

since Dgppc{0, t'p) and Dggpc{s'q, 0) vanish in our chosen coordinates and —Dpgc{0, 0) is the 
identity matrix. Due to the tensorial nature of the cross-curvature (12.11) . Co depends on 
||c||(74([/xy) and the bi-Lipschitz constants f3^ from (|4.2|) . 

Prom the above inequality and the definition of Eg^^ we deduce 

Pn = {p, e„ + ^) - {p, ^) < e\p\ + C7oe|p|2 

\q\ \q\ 

so p £ E'q as desired. □ 



6 The Monge- Ampere measure dominates the c-Monge- Ampere 
measure 

In this section we shall prove that — up to constants — the ordinary Monge- Ampere mea- 
sure \du\ dominates the c-Monge- Ampere measure \d^u\, when defined in the coordinates 
introduced in Theorem 14.31 Let us begin with a lemma which motivates our proposition 
heuristically. The conclusions of the lemma extend easily from smooth to non-smooth func- 
tions by an approximation argument combining Lemma [3TT] (b)-(c) with results of Trudinger 
and Wang [44]. However this approach would require the domains U and V to be smooth, so 
in Proposition 16.21 we prefer to construct an explicit approximation which proves the state- 
ment we need, requires no additional smoothness hypotheses, and is logically independent of 
both Lemma ETT] and [44] . 

Lemma 6.1. Assume (B0)-(B3) and let u : Uy i — > ¥i be a convex c-convex function as in 
Theorem \4-3\ If u £ C'^{U~) for some open set U'y C Uy, then \d'^u\ < 7^|(9tt| on U~, where 
7^ = ifiU^ X V) are defined as in i4-3\ )- 
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Proof. In addition to the convexity of u{q), for any y G V Theorem 14.31 asserts the convexity 
of the c-convex function q £ Uy i — > —c{q, y). Thus 

det{DlgU{q) + DggC{q,y)) < det DggU{q) 

by the concavity of S i — > det^/"(S') on symmetric non-negative definite matrices. On the 
other hand, for c-convex u G C'^{U'y)^ the measure \d'^u\ is absolutely continuous, with 
Lebesgue density given by the left hand side of (13. 7p . Thus at any q £ U^, 

<^ = det(P.V.B) + Z>;.cfeG@)) ^ 
<(^» |detD2„a(!,,G(9))| 

as desired. □ 

We now prove the proposition that we actually need subsequently. 

Proposition 6.2 (Monge- Ampere measure dominates c-Monge- Ampere measure). Assume 
(B0)-(B3), and let u : Uy i — > R be a convex c-convex function from Theorem \4-3\ Then 
{d'^ul < 7j|9n| on U'- C Uy, where -ff = jfiU- x V) are defined as in ( f^.3| j. 



Proof. It suffices to prove \d'^u\(^Br{q)) < 'y^\d'u\[Br{q)) for each ball whose closure is con- 
tained in Ug. Given such a ball, let h{q) := Ig — ^|, and let Pe{q) = e~^p{q/e) > be a smooth 
mollifier vanishing outside -Be(O) and carrying unit mass. For e > sufficiently small, we can 
define the smooth convex function 

Ue,S = {u + 5h) * Ps 



on Br{q). Since u and h are locally Lipschitz, letting R denote a bound for Lip(u) -|- Lip(/i) 
inside Br{q) yields 

\\ue,5 -{u + < eR 

for all (5 < 1. 

Claim: Fix < t < 1 and < (5 < 1. For e > sufficiently small, we claim that to each 
yo £ d^u{Btr{q)) corresponds some qe^s £ Br{q) such that {qe,&,yo) S d^Us^s- 
Indeed, for yo £ d'^u{qo) with \qo — q\ < tr, observe 

u{q) > -c{q, yo) - u"* (yo) M q£Uy, 

with equality at go- Moreover /i(g'o) < t'r- Therefore 

Ue,s{(lo) - {-c{qo,yo) -u^'iyo)) < eR + u{qo) + 6h{qo) - {-c{qo,yo) - u"* [yo)) 

= eR + 6h{qo) <eR + tr6. 

On the other hand, if g G dBr{q), then h{q) = r, and so 

Ue,s{q) - {-c{q, yo) - u""' (yo)) > -^R + u{q) + Sh{q) - {-c{q, yo) - u""* {yo)) 

> -eR + 5h{q) = -eR + r5. 
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Thus, if for fixed 6 small we choose e small enough so that 

r6 - eR> trS + eR, 

we deduce that if we lower the graph of the function —c{q, yo) — (yo) to the lowest level at 
which it intersects the graph of Ug^s, then the point of intersection must lie over Br{q). This 
proves the claim. 

Having established the claim, let E C Bj.{q) denote the (Borel) set of all q^^s which arise 
from dus^siBtriq)) in this way. Since Ue.s is smooth, the condition yo G d'^Us^siQe,5) implies 

DqUs,s(.qs,5) = -Dqc{qs,s,yo), (6.2) 

as well as 

DqqU^,s{qe,s) > - D qqC{q^ , Vo) ■ (6.3) 

By (16.21) and (Bl) we can define a smooth map Ge^siQo) throughout E using the relation 

DqUs,5{qo) = -Dgc{qo,Ge,5{qo)), 
and find that d'^Ue,5iQe,5) = {Gs,5{qE,5)} is a singleton. In this way we obtain 

\d%\{Btrm < \d'Ue,5\{E) = [ \det DgGe.5\{q)dq 

Je 

det{DggUe,5{q) + Dggc{q, Ge,5{q))) 



det D^gyCiq,Ge,5m 



dq, 



where in the last equality we used (16. 3p to deduce that DggUg^^siQ) + D gg c{q, Gs, 5 (q)) is non- 
negative definite. Hence the inequality 

det{DqqUe,s{q) + Dqqc{q,Gs,s{q))) ^ -,,n ~ i\ 
|detI.,V('?,G,,(,))| ^ detD,,n,,(,) 

holds (similarly to (|6.ip above), and so 

\d'u\{Btrm < 1~ f det{DggU,4q))dq < 7," l^ti^.^l (i?r (g))- 
Je 

Letting first e — > and then (5 — > 0, we finally deduce 



\d'u\{Btr{q)) < limsnpj^\du,,s\{Br{q)) < ^^\du\{Br{q)) . 

Here, to see the last inequality one may, for instance, use Lemma lSTTl fb) with c{x, y) = —{x, y). 
Arbitrariness of < t < 1 yields the desired result. 

□ 
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7 Alexandrov type estimates and afRne renormalization 



In this section we prove the key estimates for c-convex potential functions which will eventu- 
ally lead to the continuity and injectivity of optimal maps. Namely, we extend Alexandrov 
type estimates commonly used in the analysis of Monge- Ampere equations (thus for the cost 
c{x,y) = —{x,y)), to general non-negatively curved cost functions. This is established in 
Lemma 17.21 (plus Proposition 16. 2|) and Lemma 17.91 These estimates are used to compare 
the infimum of c-convex function on a section with the size of the section, which are the 
key ingredients in the proof of our main results; see Propositions 17.31 and 17.101 Lemma 17.91 
represents the most nontrivial and technical result we obtain in this section. 

We recall a basic lemma for convex sets due to Fritz John [24], which will play an essential 
role in the rest of the paper. 

Lemma 7.1 (John's lemma). For a compact convex set S C R", there exists an affine 
transformation L : R" R" such that Bi C L~^{S) C Bn- 

We now estimate the infimum of u in terms of the Monge- Ampere measure in a section. 
The following lemma is a standard fact for convex functions. With Lemma mi in mind, we 
state it for normalized functions u* and sections Z* . However, the estimate (17. ip is invariant 
under the affine renormalization (14.71) : according to (14.91) . it holds with or without stars. 

Lemma 7.2 (Upper bound on Dirichlet solutions to Monge- Ampere inequalities). Let u* : 

R"" ! — > RU {+00} he a convex function whose section Z* := {u* < 0} satisfies Bi C Z* C 
Bn- Assume that u* = on dZ* . Then, for all t G (0, 1), 

I ) < (l_,)n ^n^Z*) ' (^-^^ 

where tZ* denotes the dilation of Z* by a factor t with respect to the origin. 

Although the proof of this result is classical (see for instance [23]), for sake of completeness 
we prefer to give all the details. 

Proof. We can assume that u*\z* ^ 0, otherwise the estimate is trivial. It is not difficult to 
prove that 

|p*| < '"f^*"!** yp* e du*{tZ*). (7.2) 
(1 — t) 

Indeed, if q* £ tZ* and p* G du*{q*), then 

{V\q - q*) < u{q) - u{q*) = \u{q*)\ Vq G dZ*, 

and taking the supremum in the left hand side among all q G dZ* (17.21) follows. Thus, since 
^"(Z*) < C7(n), we conclude 

\du*\{tZ*)< ^^'^ 



(1-t)" ^"(Z*) 

□ 
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Combining the above lemmas, we obtain: 

Proposition 7.3. Assume (B0)-(B3) and define 7^ = 7^(Z x V) as in ll4-3\ ). Any convex 
c- convex function u : Uy 1 — > H from Theorem \4-3\ which satisfies \d'^u\ £ [A, 1/A] in a section 
of the form Z := {q G Uy\ u{q) < 0}, and u = on dZ , also satisfies 

< C(n)^|infnp. (7.3) 
A z 

Proof. First use the affine map L as given in Lemma [7. II to renormahze u into u* using (14. 7p . 
This does not change the bound \du*\ G [A, 1/A], but allows us to apply Lemma [72] with 
t = 1/2. Its conclusion (|7.ip has been expressed in a form which holds with or without the 
stars, in view of (|4.9I) . Proposition 16.21 now yields the desired inequality (|7.3p . □ 

7.1 c-cones over convex sets 

We now progress toward the Alexandrov type estimate in Lemma \7M In this subsection we 
construct and study the c-cone associated to the section of a c-convex function. This c-cone 
plays an essential role in our proof of Lemma 17.91 

Definition 7.4 (c-cone). Assitme (B0)-(B2) and (A3w), and let u : Uy 1 — > R be the 

c-convex function with convex level sets from Theorem \4.3[ Let Z denote the section {u < 0}, 
fix q £ intZ, and assume u = on dZ. The c-cone h'^ : Uy 1 — > R generated by q and Z 
with height —u{q) > is given by 

h^iq) := sup{-c(g, y) + c{q, y) + u{q) \ -c{q, y) + c{q, y) + u{q) < on dZ}. (7.4) 
ydV 

Notice the c-cone hf^ depends only on the convex set Z C Uy, q £ int Z, and the value 
u{q), but is otherwise independent of u. Recalling that c{q,y) = on ?7y, we record several 
key properties of the c-cone: 

Lemma 7.5 (Basic properties of c-cones). Adopting the notation and hypotheses of Definition 
\7.4\ let h'^ : Uq 1 — > R be the c-cone generated by q and Z with height —u{q). Then 

(a) h'^ has convex level sets; it is a convex function if (B3) holds; 

(b) h^{q) > h^{q) = u{q) for all q G Z; 

(c) h^ = on dZ; 

(d) d^h^iq) C d^u{Z). 

Proof. Property (a) is a consequence of the level-set convexity of q 1 — > —c{q, y) proved in 
Theorem 14.31 or its convexity assuming (B3). Moreover, since —c{q,y)+c{q,y)+u{q) = u{q) 
for all q £ Uy, (b) follows. For each pair z £ dZ and yz G d'^u{z), consider the supporting 
mountain mz{q) = —c{q,yz) + c{z,yz), i.e. mz{z) = = u{z) and niz < u. Consider the 
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c-segment a{t) connecting cj(0) = y and cj(1) = yz in V with respect to z. Since —c{q, y) = 0, 
by continuity there exists some t E]0, 1] for which rhz{q) := —c{q,a{t)) + c{z,a{t)) satisfies 
fhz{q) = u{q). Prom Loeper's maximum principle (Theorem 13.31 above), we have 

rhz < inax[mz,—c{-,y)] = m.ax[mz,0], 

and therefore, from < tt, 

rhz < on Z. 

By the construction, rhz is of the form 

-c{-,y) +c{q,y) + u{q), 

and vanishes at z. This proves (c). Finally (d) follows from (c) and the fact that h'^{q) = u{q). 
Indeed, it suffices to move down the supporting mountain of at q until the last moment at 
which it touches the graph of u on Z from below. The conclusion then follows from Loeper's 
local to global principle, Corollarv 13.41 above. □ 

The following estimate shows that the Monge- Ampere measure, and the relative location 
of the vertex within the section which generates it, control the height of any well-localized 
c-cone. Afficionados of the Monge- Ampere theory may be less surprised by this estimate once 
it is recognized that the localization in coordinates ensures the cost is approximately affine, 
at least in one of its two variables. Still, it is vital that the approximation be controlled! 
Together with Lemma iTTST d). this proposition plays a key role in the proof of our Alexandrov 
type estimate (Lemma 17. 9p . 

Proposition 7.6 (Lower bound on the Monge- Ampere measure of a small c-cone). Assume 
(B0)-(B3) and define c £ C^{Uy x V) as in Definition \4.1\ Let Z C Uy be a closed convex 
set and the c-cone generated by q £ mi Z of height —h^{q) > over Z. Let n"*",!!" 
be two parallel hyperplanes contained in T^V \ Z and touching dZ from two opposite sides. 
Then there exists £c > small, depending only on the cost (and given by Lemma m\) . and a 
constant C{n) > depending only on dimension, such that if diam{Z) < ec/C{n) then 

mr < c(n) "^^"^^"^(^''"/)'^"^(^''"'» ia/.^i({g})^"(z), (7.5) 

where £n+ denotes the maximal length among all the segments obtained by intersecting Z with 
a line orthogonal to . 

To prove this, we first observe a basic estimate on the cost function c. 

Lemma 7.7. Assume (B0)-(B2). For c G C^{Uy x V) from DefinitionlJJ] and each y GV 
and q,q & Uy, 

I - L)gc{q, y) + Dqc{q, y)\ < —\q - q\ \Dqc{q, y)\ (7.6) 
where Ec is given by = '^{fit)^{Pc)^\\-^'xxy'^\\L°°{UxV) notation (|4.2p . 
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Proof. For fixed q £ Uy introduce the c-exponential coordinates p{y) = —Dqc{q,y). The 
bi-Lipschitz constants (14. 2p of this coordinate change are estimated by f3f < P'^ I3~ as in 
Corollary [131 Thus 

dist(y, y) <(3^\- Dqc{q, y) + Dqc{q, y)\ 
= (3tf3~\Dqc{q,y)\. 

where c{q,y) = from Definition 14.11 has been used. Similarly, noting the convexity (B2) of 
Vii:=p{V), 

I - Dqc{q, y) + Dqc{q, y)\ = \- Dgc{q, y) + Dqc{q, y) + Dqc{q, y) - Dqc{q, y)\ 

< \\DlqDpc\\^^(^^_^y_^\q - q\\p{y) - p{y)\ 

< \\DlqDyc\\Lc^^u.^v){l3c l3t?\Q - q\ dist(y, y) 

The result follows since \D^qDyc\ < {{/3-f + P+{f3-)'-^))\Dl^Dyc\ < 2(3+ {(3-f\Dl^Dyc\ . 
(The last inequality follows from (3~ > 1.) □ 

Proof of Proposition 1 7. 61 We fix g G Z. Let IT*, i = 1, (with equal either 11"^ or 

n~) be hyperplanes contained in R" \ Z, touching 5Z, and such that {11"'", 11^, . . . ,11'^} are 
all mutually orthogonal (so that also {n~,n^, ... ,11"} are mutually orthogonal). Moreover 
we choose {11^, . . . ,11"} in such a way that, if tt^{Z) denotes the projection of Z on 11^ and 
"'~^(7r^(Z)) denotes its (n — l)-dimensional Hausdorff measure, then 

n 

C(n)jr""-i(7ri(Z)) > []dist(g,ff), (7.7) 

i=2 

for some universal constant C{n). Indeed, as vr^(Z) is convex, by Lemma \T7\\ we can find an 
ellipsoid E such that E C vr^(Z) C (n — 1)E^ and for instance we can choose {n^, . . . jLt"} 
among the hyperplanes orthogonal to the axes of the ellipsoid (for each axis we have two 
possible hyperplanes, and we can always choose the furthest one so that (|7.7p holds). 

Each hyperplane 11* touches Z from outside, say at q^ E T^V . Let pi G TyV be the 
outward (from Z) unit vector at orthogonal to 11*. Then sipi G dh^{q^) for some Sj > 0, 
and by Corollary 13.41 there exists yi G d^h'^[q^) such that 

-Dqc{q\yi) = SiPi. 

Define yi{t) as 

-Dqc{q\yi{t)) = tsipi, 

i.e. yi{t) is the c-segment from y to yi with respect to g*. As in the proof of Lemma [731 (c), 
the intermediate value theorem yields < < 1 such that 

-c{-,y,{ti)) + c{q,yi{ti)) + h^{q) < on Z 
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Figure 2: The dotted line represents tiie grapii of — c(-, y^) + c{q, yi) + h'^{q), while the dashed 

one represents the graph of mi{ti) := —c{-,yi{ti)) + c{q,yi{ti)) + h'^{q). The idea is that, whenever 
we have a supporting function for /i^ at a point € dZ , we can let y vary continuously along the 
c-segment from y to yi with respect to g\ to obtain a supporting function rriiiti) which touches h'^ 
also at y. 



with equality at g*. Thus, by the definition of h'^, yiiU) G d'^h'^{q) n 

-Dqc{q,yi{ti)) G and tiSip, = -Dqc{q\yi{ti)) G dh''{q'). 

Therefore by the convexity of /i"^ shown in Lemma [73] (a). the aflfine function with slope 
-Dqc{q\yi{ti)) and with P*(n*) = satisfies P\q) < h%q). This shows 



dist(g, m)' 

Also, by (TTTHIl 



I - Dqc{q,yi{ti)) + Dqc{q\yi{ti))\ < — |g - | - L'gc(g, 

< — diamZ| - Dqc{q,yi{ti)\. 

Therefore if diamZ < 5n£c with 5„ > small, each vector —Dqc{q,yi{ti)) is close to 
-Dqc{q\yi{ti)), say 

I - Dqc{q,yi{ti)) + i:>gc(g*,yi(tj))| < - Dqc{q,yi{ti))\. 

Since the vectors {—Dqc{q^ ,yi{ti))} are mutually orthogonal, the above estimate implies 
that for 6n small enough the convex hull of {— Dqc{q, yi(ti))} has measure of order HlLi I ~ 
Dqc{q^ ,yi{ti))\. Thus, by the lower bound (17. Sh and the convexity of dh'^{q), we obtain 

^""{dh^iq)) > C{n) ' 



niLidist(g,m)- 
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Figure 3: The volume of any convex set always controls the product (measure of one sHce) • (measure 
of the projection orthogonal to the slice). 



Since 11 was either 11+ or 11 , we have proved that 

n 

\h^{qT <C{n) min{dist (q', n+ ) , dist (g, ) } J]^ dist (g, Jl')\dh^\{{q]) . 

i=2 

To conclude the proof, we apply Lemma 17.81 below with Z' given by the segment obtained 
intersecting Z with a line orthogonal to 11"'". Combining that lemma with (|7.7I1 . we obtain 

n 

C{n)\Z\ > eu+Yldistiq,^), 

i=2 

and last two inequalities prove the proposition (taking C{n) > l/6n larger if necessary). □ 

Lemma 7.8 (Estimating a convex volume using one slice and an orthogonal projection). Let 
Z be a convex set in R" = R" x R** . Let 7r',7r" denote the projections to the components 
R" , R" , respectively. Let Z' be a slice orthogonal to the second component, that is 

Z' = {tt")~^{x") n Z for some x" G t^'\Z). 

Then there exists a constant C{n), depending only on n = n' + n" , such that 

C{n)^''{Z) > ^"'(Z')jr""(7r"(^)), 

where M"^ denotes the d-dimensional Hausdorff measure. 
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Proof. Let L : R" — > R" be an afRne map with determinant 1 given by Lemma 17.11 such 
that Bj- C L{tt"{Z)) C Bn"r for some r > 0. Then, if we extend L to the whole R" as 
L{x',x") = {x',Lx"), we have ^"(L(Z)) = jr"'(L(Z')) = and 

^""(7r"(L(Z))) = ^""(L(7r"(Z))) = 

Hence, we can assume from the beginning that Br C vr"(Z) C i?n"r- Let us now consider 
the point x" , and we fix an orthonormal basis {ei, . . . ,e„"} in R" such that x" = cei for 
some c < 0. Since {rei, . . . ,re„"} C 7r"(Z), there exist points {xi, . . . ,Xn"} C Z such that 
7r"(xj) = rej. Let C denote the convex hull of Z' with xi, and let V denote the (n' + 1)- 
dimensional strip obtained taking the convex hull of R" x {x"} with xi. Observe that 
C cV, and so 

jr"'+i(C") = ^^^dist(xi,R"' X > -^^"'(Z'). (7.9) 

We now remark that, since '^"{xi) = rei and _L 1^' for i = 2, . . . , n", we have dist(xj, V') = r 
for alH = 2, . . . ,n". Moreover, if yi G denotes the closest point to Xj, then the segments 
joining Xi to yi parallels e^, hence these segments are all mutually orthogonal, and they are 
all orthogonal to V too. From this fact it is easy to see that, if we define the convex hull 

C := co{x2, . . . ,x„/',C'), 

then, since \xi — yi\ = r for i = 2, . . . , n" , by (17. 9p and the inclusion vr"(Z) C Bn"r C R"" we 
get 

This concludes the proof, as C C □ 
7.2 An Alexandrov type estimate 

The next Alexandrov type lemma holds for localized sections Z of c-convex functions. 

Lemma 7.9 (Alexandrov type estimate and lower barrier). Assume (B0)-(B3) and let 

u : Uy I — > H be a convex c-convex function from Theorem \4-3\ Let Z denote the section 
{u < 0}, assume u = on dZ, and fix q £ intZ. Let 11^,11" be two parallel hyperplanes 
contained in R" \ Z and touching dZ from two opposite sides. Then there exists e'f.{n) > 
small, depending only on dimension and the cost function (with £'^{n) = £c/C{n) given by 
Proposition 1 7. d\) such that if d.\axn.{Z) < e'^{n) then 

mr < CinUiZ X V) "^"^idist(g, n+), dist(g, H-)} ^z)^n^z), 

where denotes the maximal length among all the segments obtained by intersecting Z with 
a line orthogonal to 11+ , and 7^ = j^{Z x V) is defined as in |^.3| ). 
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Proof. Fix q G Z. Observe that n = on dZ and consider the c-cone h'^ generated by q and 
Z of height —h^{q) = —u{q) as in (17. 4p . From Lemma [731 (d) we have 

\d'h'\m)<\d'u\{z), 

and from Loeper's local to global principle, Corollary 13.41 above. 

dh%q) = -D,c{q,d'h%q)). 

Therefore 

\dh'\m) < \\detDlc\\co^{g}^v)\d'h'^\m)- 
The lower bound on comes from (|7.5|) . This finishes the proof. □ 



7.3 Estimating solutions to the c-Monge- Ampere inequality \d'^u\ G [A,l/A] 
Combining the results of Proposition 17.31 and Lemma ED yields: 

Proposition 7.10 (Bounding local Dirichlet solutions to c-Monge- Ampere inequalities). 
Assume (B0)-(B3) and let u : Uy i — > R &e a convex c-convex function from Theorem \4-3[ 
There exists e'^{n) > small, depending only on dimension and the cost function (and given 
by Lemma \7.9\} . and constants C{n), Ci{n) > 0, i = 1,2, depending only on dimension, such 
that the following holds: Letting Z denote the section {u < 0}, assume \d^u\ G [A, 1/A] in Z 
and u = on dZ . Let 11+ ^ 11" be parallel hyperplanes contained in T^V\Z and supporting 
Z from two opposite sides. //diam(Z) < e'^{n) then 



and 



l%)r <c(,,)2£"-Mdist(,,n+),dist(,,n-)} ^^^.^^^^^ ^^^^^^ 



^"(Z)2 - ^ ^ A iu+ 
where £yi+ denotes the maximal length among all the segments obtained by intersecting Z with 
a line orthogonal to 11"'', and 7^ = ^^{Z x V) is defined as in (J^- 

Proof. Equation (|7.1ip follows from Lemma [7^ and the assumption \d'^u\ < 1/A. Now, by 
Lemma 17.11 we deduce that there exists an ellipsoid E such that E C Z C nE, where nE 
denotes the dilation of by a factor n with respect to its barycenter q. Taking 11+ and IT" 
orthogonal to one of the longest axes of E and q = q in (17.111) yields 

imr < c{n\^\j£-{zf. 

On the other hand, convexity of u along the segment which crosses Z and passes through 
both q and the point q where inf^ u is attained implies 

I inf u\ < n\u{q)\, 

since the barycenter of E divides the segment into a ratio at most n : 1. Combining these 
two estimates with (17. 3p we obtain (l7.1Qp . to complete the proof. □ 
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Remark 7.11 (Stability of bounds under afRne renormalization). Noting 7g < 7(^7(7 from 
Corollary 14.41 we observe that the estimate (I7.1QP is stable under affine renormalization: let 
L be an afRne transformation, and recall the renormalization 

u*{q) := |detLr2/"i2(Lg). 

Then L~^{Z) is a section for u* and 

I inf u*[' = \deiLr\miu\'' ^\deiL\~^^'\Zf = ^'\L'\Z)f. 

On the other hand, estimate (|7.11l) is not stable under affine renormalization (a line orthog- 
onal to 11"'" is not an affinely invariant concept). For this reason, both in the proof of the 
c-strict convexity (Section [8|) and in the proof of differentiability u G (Section (H) we ap- 
ply our Alexandrov estimates directly to the original sections, without renormalizing them. 
Using this strategy, our estimates turn out to be strong enough to adapt to our situation the 
strict convexity and interior continuity theory of Caffarelli [5j [8j. We perform this in the 
remainder of the manuscript. 



8 The contact set is either a single point or crosses the domain 

In this section and the final one, we complete the crucial step of proving the strict c-convexity 
of the c-convex potentials u : U t — > R arising in optimal transport, meaning d'^u{x) should be 
disjoint from d'^u{x) whenever i € U'^ are distinct. This is accomplished in Theorem 19.11 
In this section, we show that, if the contact set does not consist of a single point, then it 
extends to the boundary of U . Our method relies on the non-negative cross-curvature (B3) 
of the cost c. 

Prom now on we adopt the following notation: a ~ 6 means that there exist two positive 
constants Ci and C2, depending on n and "^t^c only, such that Cia < b < C2a. Analo- 
gously we will say that a < b (resp. a > 6) if there exists a positive constant C, depending 
on n and ^^Ic only, such that a < Cb (resp. Ca > b). 

Recall that a point x of a convex set S C R" is exposed if there is a hyperplane supporting 
S exclusively at x. Although the contact set S := d'^* u'^* (y) may not be convex, it appears 
convex from y by Corollary 13.41 meaning its image q{S) C Uy in the coordinates (|4.4p is 
convex. The following theorem shows this convex set is either a singleton, or contains a seg- 
ment which stretches across the domain. We prove it by showing the solution geometry near 
certain exposed points of q{S) inside Uy would be inconsistent with the bounds established 
in the previous section. 

Theorem 8.1 (The contact set is either a single point or crosses the domain). Assume 
(B0)-(B3) and let u be a c-convex solution of (|4.ip with C U open. Fix x ^ and 
y G d^u{x), and define the contact set S := {x £ U \ u{x) = u{x) — c(x, y) + c{x, y)}. Assume 
that S 7^ {x}, i.e. it is not a singleton. Then S intersects dU . 
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Figure 4: If the contact set Sy has an exposed point qo, we can cut two portions of Sy with two 
hyperplanes orthogonal to q — (7°. The diameter of —Dyc{Ko,y) needs to be sufficiently small to 
apply the Alexandrov estimate Lemma FflOl while —Dyc{KQ,y) has to intersect is a set of positive 
measure to make use of Lemma [731 in the case g° is not an interior point of spt \d^u\. 

Proof. As in Definition 14.11 we transform {x,u) 1 — > {q,u) with respect to y, i.e. we consider 
the transformation q & Uy 1 — > x{q) G U, defined on Uy := —Dyc{U,y) C T^V by the 
relation 

-Dyc{x{q),y) = q, 

and the modified cost function c{q,y) := c{x{q),y) — c{x{q),y) on Uy x V, for which the 
c-convex potential function q £ Uy 1 — > u{q) := u{x{q)) + c{x{q),y) is convex. We observe 
that c{q, y) = for all q, and moreover the set S = d'^ u'^ (y) appears convex from y, meaning 
Sy := —Dyc{S,y) C Uy is convex, by the Corollary 13.41 to Loeper's maximum principle. 

Our proof is reminiscent of Caffarelli's for the cost c{q, y) = —{q, y) [HI Lemma 3]. Observe 
q := —Dyc{x,y) lies in the interior of the set U^ := —Dyc{U^,y) where \d'^u\ S [A/7^,7~/A], 
according to Corollary 14.41 Choose the point q^ £ Sy C Uy furthest from q; it is an exposed 
point of Sy. We claim either q^ = q or q^ dUy. To derive a contradiction, suppose the 
preceding claim fails, meaning q^ £ Uy\ {q}. 

For a suitable choice of Cartesian coordinates on V we may, without loss of generality, 
take q^ — q parallel to the positive y^ axis. Denote by the associated orthogonal basis for 
TyV, and set 6° := (g°,ei) and b := {q,ei), so the halfspace qi = {q,ei) > 6° of T^V ~ R" 
intersects 5^;^ only at g". Use the fact that is an exposed point of Sy to cut a corner Kq 
off the contact set 5 by choosing s > small enough that 6 = (1 — s)b^ + sb satisfies: 

(i) —Dyc{KQ,y) := Sy f] {q £ Uy \ qi > b} is a compact convex set in the interior of Uy; 

(ii) d[am{—Dyc{Ko,y)) < ec/2, where is from Lemma [7?9l 
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I 

Figure 5: We cut the graph of u with the two functions and ml to obtain two sets Kg « Kg and 
Kl w Kq inside which we can apply our Alexandrov estimates to get a contradiction (both Lemma 
17.21 and Lemma FflOl to K^, but only Lemma [731 to K^: ). The idea is that the value of u — ml at xq is 
comparable to its minimum inside i^e, but this is forbidden by our Alexandrov estimates since xq is 
too close to the boundary of Kq. However, to make the argument work we need also to take advantage 
of the section K^:, in order to "capture" some positive mass of the c-Monge- Ampere measure. 

Defining := (1 — s)q^ + sq, := x{q'^) the corresponding c-segment with respect to 
y, and q = q^, note that Sy n {qi = b} contains q, and Kq contains x := and 
Since the corner Kq needs not intersect the support of \d'^u\ (especially, when q^ is not 
an interior point of spt|(9'^u|), we shall need to cut a larger corner Kq as well, defined by 
—Dyc{KQ,y) := Sy f^ {q ^ U y \ qi > 6}, which intersects U'^ at x. By tilting the supporting 
function slightly, we shall now define sections i^^ C Kl of u whose interiors include the 
extreme point x^ and whose boundaries pass through x and x respectively, but which converge 
to Kq and Kq respectively as e — > 0. 
Indeed, set := y + eei and observe 

ml{x) := -c{x, ye) + c(x, y) + c(x^ y^) - c(x^ y) 
= e{-Dyc{x, y) + Dyc{x'',y), ei) + o(e) 

= e{{-Dyc{x, y), ei) - (1 - - sb) + o{e) (8.1) 

Taking s G {s, 1} in this formula and e > shows the sections defined by 

Ke ■■= {x I u{x) < u{x) - c{x, ye) + c(x, ye)}, 
Kl := {x I u{x) < u{x) - c{x, ye) + c(x, ye)}, 

both include a neighbourhood of xq but converge to Kq and Kq respectively as e — > 0. 

We remark that there exist a priori no coordinates in which all set Ke are convex. However 
for each fixed e > 0, we can change coordinates so that both Ke and K^ become convex: use 
ye to make the transformations 

q := -Dyc{xe{q),ye), 

Ceiq,y) ■■= c{Xe{q),y) - c{Xe{q),ye), 
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so that the functions 



Usiq) ■■= u{xe{q)) + c{xeiq), Ve) - u{x) - c(x, y^), 
ul{q) := u{xe{q)) + c{xe{q), Ve) - u{x) - c{x, Ve). 

are convex on Uy^ := Dyc{U,ye). Observe that, in these coordinates, and become 
convex: 

:= -Dyc{K„y,) = {q e Uy^ \ uM < 0}, 
kl := -Dyc{Kly,) = {g G Uy^ \ ul{q) < 0}, 

and either C k^ or k^ C since Ue(g) — ul{q) = const. For e > small, the inclusion 
must be the first of the two since the limits satisfy ko C -ftTg and q £ k^ \ k^. 

In the new coordinates, our original point x G C/^, the exposed point x^ , and the c-convex 
combination x with respect to y, correspond to 

Qe ■■= -Dyc{x,ye), ql := -Dyc{x^,ye), qs := -Dyc{x,ys). 

Thanks to (ii), we have diam(i^e) < for e sufficiently small, so that the estimate of Lemma 
17.91 applies. In these coordinates (chosen for each e) we consider the parallel hyperplanes 
n+ / nj which support kgr C Uy^ from opposite sides and which are orthogonal to the 
segment joining q^ with q^. Since lime^o'Ze ~ Qe = ~ Q paralles the ei axis the limiting 
hyperplanes 11^ = lime_>on^ must coincide with = {q £ TyV \ qi = and IIq = {q £ 
TyV \qi = h}. Thus g° G n+ and 

Observing that \q^ — qe\ is shorter than segment obtained intersecting k^ with the line 
orthogonal to and passing through q^ G int k^ , Lemma 17.91 combines with C and 
\d^^Ue\{K^) < A-f-^'^iK,) from gl|) and Corollary |131 to yield 

-0 as.^0. (8.2) 



A7,-^"(Ki)2 

On the other hand, x e S imphes Ue{q^) = —m^{x^), and x G 5 implies uliqs) = —ml{x^) 
similarly. Thus (|8.ip yields 

uM) e{W-b) + o{e) ^ ' 

Our contradiction with (|8.2|) - (IO]l will be estabhshed by bounding the ratio (q^)]"^ / J^'^ (Kl)"^ 
away from zero. 
Recall that 

}p = {-Dyc{x^ ,y),ei) = m.&y.{qi\ q £ -Dyc{KQ,y)}) > b 
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and u{x) — u{x) > — c(x, y) + c(x, y) with equality at x^. Prom the convergence of to 
and the asymptotic behaviour (|8.ip of ml{x) we get 

uliq^) -u{x^) -c{x^,ye) + u{x) + c{x,ye) 



inf^i^i supgg^i[-u(x(g)) -c{x{q),ys) + u{x) +c{x,y^)] 
-c(x°, ye) + c(x, ye) + c(a;°, y) - c(x, y) 



> 



> 



snPxGK]: [-cix, ye) + c{x, ye) + c(x, y) - c{x, ] 
e{{-Dyc{x^,y), ei) -h) + o{e) 

e(max{gi | q G -Dyc{K^,y)} -b) + o{e) 
1 



(8.4) 



> 
- 2 

for e sufficiently small. This shows u^iqe) is close to the minimum value of ■u^. We would like 
to appeal to Lemma 17.21 to conclude the proof, but are unable to do so since we only have 
bounds \d'^u\ G [A, 1/A] on the potentially small intersection of with Ke- However, this 
intersection occupies a stable fraction of as e — > 0, which we shall prove as in [H Lemma 
3]. 

Since converges to Kq for sufficiently small e, observe that is uniformly bounded. 
Therefore the affine transformation {Ll)~^ that sends Ke to Bi C Ke'* C Bn as in Lemma 
17. H is an expansion, i.e. \{Ll)~^q— {Ll)~^q'\ > Co\q — q'\, with a constant Cq > independent 
of £. Since x is an interior point of U^, B^^- ^i^^{x) C for sufficiently small (5 > 0, hence 

B2S/CoiQe) C U^^ with (5- from gj). Defining U^^ := -Dyc{U^ ,ye), we have 

Uy)* := {Ll)-\UI) D B2s{q*e). 
Reducing 6 if necessary to ensure 6 < 1, define (to apply Lemma [72] later) 

Kl;; := (1 - 5)K^'*. 

(As in Lemma [721 (1 ~ S)Ke'* denotes the dilation of Ke'* of a factor (1 — 6) with respect 
to the origin.) Since Ke'* is convex, it contains the convex hull of Bi U {^g}, and so 

r>l,*\ 



^^{B^smnK '*s)>C5^. 



for some constant C = C{n) > depending on dimension only. Letting K^ ^ := L].{K^^ 
this implies 

^'\Ul n Kl^) > C\ detL^|(5" ~ ^"(K^i)5". 

Recalling that > and < denote inequalities which hold up to multiplicative constants de 
pending on n, A and ^tlt /^j Proposition 16.21 combines with this estimate to yield 

\du\\{klA > \d'^ul\{K}s) > ^""{kDS^ 
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where (14.11) and Corollary 14.41 have been used. Finally, since the conclusion of Lemma 17.21 
holds with or without stars in light of (|4.7p - (l4.9p . taking t = (1 — J) in (|7.ip yields 



inf^i M^r 



Since (5 > is independent of £ this contradicts (|8.2p - (|8.3p to complete the proof. □ 

Remark 8.2. As can be easily seen from the proof, one can actually show that \i = U 
and S is not a singleton, then Sy has no exposed points in the interior of Uy. Indeed, if by 
contradiction there exists an exposed point of Sy belonging to the interior of Uy, we can 
choose a point q ^ Sy \n the interior of Uy = U^ such that the segment q^ — q\s orthogonal 
to a hyperplane supporting 5^; at q^ . Then it can immediately checked that the above proof 
(which could even be simplified in this particular case) shows that such a point q^ cannot 
exist. 



9 Continuity and injectivity of optimal maps 

The first theorem below combines results of Sections [5] and [8] to deduce strict c-convexity of 
the c-potential for an optimal map, if its target is strongly c-convex. This strict c-convexity 
— which is equivalent to injectivity of the map — will then be combined with an adapta- 
tion of Caffarelli's argument [H Corollary 1] to obtain interior continuity of the map — or 
equivalently C^-regularity of its c-potential function — for non-negatively cross-curved costs, 
yielding the concluding theorem of the paper. 

Theorem 9.1 (Injectivity of optimal maps to a strongly c-convex target). Let c satisfy (BO)- 
(B3) and (B2u). If u is a c-convex solution of (|4.ip on U^ <Z U open, then u is strictly 
c-convex on U'^ , meaning d'^u{x) and d'^u{x) are disjoint whenever x,x € U^ are distinct. 

Proof. Suppose by contradiction that y G d^u{x) n d'^u{x) for two distinct points x,x € U'^, 
and set S = (y). According to Theorem 18.11 the set S intersects the boundary of U 

at a point x G dU Ci d'^*u^*{y). Since ()4.ip asserts A < \d'^u\ on and \d'^u\ < A on U, 
Theorem IS.lT a) yields y (since x,x £ C/^), and hence x £ U hy Theorem 15.11 (b). This 
contradicts x € dU and proves the theorem. □ 

Theorem 9.2 (Continuity of optimal maps to strongly c-convex targets). Let c satisfy (BO)- 
(B3) and (B2u). Ifu is a c-convex solution of (|4.ip on U'^ C U open, then u is continuously 
differentiahle inside U^ . 

Proof. Recalling that c-convexity implies semiconvexity, all we need to show is that the c- 
subdifferential d^u(x) of u at every point 5; G C/a is a singleton. 

Assume by contradiction that is not. As d^u{x) is compact, one can find a point yo in 
the set d^u{x) such that —Dxc{x,yo) S du{x) is an exposed point of the compact convex 
set du{x). Similarly to Definition 14.11 we transform (x,n) i — > {q,u) with respect to yo, i.e. 
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Figure 6: w e du{0) and the hyperplane orthogonal to v is supporting i9u(0) at 0. 



V 




I H 

~e/o(l) ~e 

Figure 7: Since the hyperplane orthogonal to v is supporting du{0) at 0, we have u{—tv) = o(t) for 
t > 0. Moreover, u grows at least linearly in the direction of v. 

we consider the transformation q G Uyg i — > x{q) G U , defined on Uy^ = —Dyc{U,yo) + 
Dyc{x,yQ) C Ty^V by the relation 

-Dyc{x{q),yo) + Z:>yc(x, yo) = q, 

and the modified cost function c{q,y) := c{x{q),y) — c{x{q),yo) on Uy^^ x V, for which the 
c-convex potential function q £ Uyf^ i — > u{q) := u{x{q)) — u{x) + c{x{q),yQ) — c{x,yQ) is 
convex. We observe that c{q, yo) = for all q, the point x is sent to 0, u > u(0) = 0, and 
u is strictly convex thanks to Theorem I9.1[ Moreover, since —Dxc{x, yo) G du{x) was an 
exposed point of du{x), = — Dgc(0,yo) is an exposed point of 5£t(0). Hence, we can find a 
vector V G du{Qi) \ {0} such that the hyperplane orthogonal to f is a supporting hyperplane 
for 9n(0) at 0. Thanks to the convexity of u, this implies that 

u{-tv) = o{t) for t > 0, u{q) > {v, q) + n(0) for all q £ Uy^,. (9.1) 

Let us now consider the section := {u < e}. Since u{0) = 0, tt > and u is strictly 
convex, Ks — > {0} as e — > 0. Thus by (|9.1I1 it is easily seen that for e sufficiently small the 
following hold: 

Ke C{q\ {q,v) < e}, -a{e)v G K^, 
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where a(e) > is a positive constant depending on e and such that a{e)/e +00 as e — > 0. 
Since is the minimum point of u, this immediately imphes that one between our Alexandrov 
estimates (|7.10|) or (17. lip must be violated by u inside for e sufficiently small, which is 
the desired contradiction. □ 
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